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Abstract. We survey rigorous, formal, and numerical results on the formation of 
point-like singularities (or blow-up) for a wide range of evolution equations. We use a 
similarity transformation of the original equation with respect to the blow-up point, 
such that self-similar behaviour is mapped to the fixed point of a dynamical system. 
We point out that analysing the dynamics close to the fixed point is a useful way 
of characterising the singularity, in that the dynamics frequently reduces to very few 
dimensions. As far as we are aware, examples from the literature either correspond 
to stable fixed points, low-dimensional centre-manifold dynamics, limit cycles, or 
travelling waves. For each "class" of singularity, we give detailed examples. 
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1. Introduction 

Non-linear partial differential equations (PDE's) are distinguished by the fact that, 
starting from smooth initial data, they can develop a singularity in finite time [DElEli]. 
m Very often, such a singularity corresponds to a physical event, such as the solution 
(e.g. a physical flow field) changing topology, and/or the emergence of a new (singular) 
structure, such as a tip, cusp, sheet, or jet. On the other hand, a singularity can also 
imply that some essential physics is missing from the equation in question, which should 
thus be supplemented with additional terms. (Even in the latter case, the singularity 
may still be indicative of a real physical event). 

Consider for example the physical case shown in Fig. [H which we will treat in 
section H] below. Shown is a snapshot of one viscous fluid dripping into another fluid, 
close to the point where a drop of the inner fluid pinches off. This process is driven 
by surface tension, which tries to minimise the surface area between the two fluids. At 
a particular point xo,to in space and time, the local radius h{x,t) of the fluid neck 
goes to zero; this point is a singularity of the underlying equation of motion. Since 
the drop breaks into two pieces, there is no way the problem can be continued without 
generalising the formulation to one that includes topological changes. However, in this 
review we adopt a broader view of what constitutes a singularity. We consider it as such 
whenever there is a loss of regularity, which implies that there is a length scale which 
goes to zero. This is the situation under which one expects self-similar behaviour, which 
is our guiding principle. 

A fascinating aspect of the study of singularities is that they describe a great 
variety of phenomena which appear in the natural sciences and beyond [3J. Some 
examples of such singular events occur in free-surface flows [6], turbulence and Euler 
dynamics (singularities of vortex tubes [71 [8] and sheets [9]), elasticity [10], Bose-Einstein 
condensates [H], non-linear wave physics [12], bacterial growth [T3l [H], black- hole 
cosmology [15], |16] , and financial markets [17] . 

In this paper we consider evolution equations 



where F[h] represents some (nonlinear) differential or integral operator. We will also 
discuss cases where h is a vector, and thus (11. ip is a system of equations. Furthermore, 
the spatial variable x may also have several dimensions, and thus potentially different 
scaling in different coordinate directions. We will cite some examples below, but few 
of the higher- dimensional cases have so far been analysed in detail. For the purpose 
of the following discussion, let us suppose that both x and h are scalar quantities, and 
that the singularity occurs at a single point in space and time xo,to- If t' = to — t and 
x' = X — Xq, we are looking for local solutions of (11.11) which have the structure 



ht = F[h], 






I Of course, there are also many examples of nonlinear PDE's for which global existence can be 
established! 
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Figure 1. A drop of Glycerin dripping through Polydimethylsiloxane near pinch-ofF 
[5]. The nozzle diameter is 0.48 cm, the viscosity ratio is A = 0.95. 

with appropriately chosen values of the exponents a,f3. Note that later the prime is 
also used to indicate a derivative. However, this will always be with respect to a spatial 
variable like x, z, or the similarity variable ^, hence confusion should not arise. 

Giga and Kohn [18], [19] proposed to introduce self-similar variables r = — In(t') 
and C, = x' /t'^ to study the asymptotics of blow up. Namely, putting 

h{x,t)=t"'H{^,T), (1.3) 

(11. ip is turned into the "dynamical system" 

Hr = G[H] = aH - (3^H^ + F[H]. (1.4) 

By virtue of (II. 4p . solutions to the original PDE (11.10 for given initial data can be viewed 
as orbits in some infinite dimensional phase phase, for instance, L^. To understand the 
blow-up of (II. ip . Giga and Kohn proposed to study the long-time behaviour of the 
dynamical system (11.40 . Thus in particular, one is interested in the attractors of (11.40 
(tu-limit sets in the notation which is customary in the context of partial differential 
equations, see [20] and references therein). If (II. 2p is indeed a solution of (II. ip . the right 
hand side of (II. 4p is independent of r, and self-similar solutions of the form (11.21) are 
fixed points of (II. 4p . which we will denote by H{^). By studying the dynamics close to 
the fixed point, we find that the dynamical system (11.41) frequently reduces to very few 
dimensions. Thus on one hand one obtains detailed information on the behaviour of the 
original problem (II. ip near blowup. On the other hand, one also gains a fruitful means 
of classifying, or at least characterising singularities. 

The most basic linear stability analysis of this self-similar solution consists in 
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linearising around the fixed point according to 



(1.5) 



which gives 



Pr = CP, 



(1.6) 



where C = C{H) depends on the fixed point solution H. To solve (11.61) . we write P as 
a superposition of eigenfunctions Pj of the operator C: 



In the cases we know, the spectrum turns out to be discrete. For evolution 
PDE's involving second order elliptic differential operators, such as semilinear parabolic 
equations, mean curvature or Ricci flows, the discreteness of the spectrum of the 
linearisation about the flxed point is a direct consequence of Sturm-Liouville theory 
[211 [22]. This theory establishes that, under quite general conditions on the coefficients 
of a second order linear differential operator and the boundary conditions, its spectrum is 
discrete and the corresponding eigenfunctions form a complete set in a suitably weighed 
space. Some explicit examples are presented in subsection 13.1.11 For general linear 
operators such a theory is not available, and one has to study the spectrum case by 
case. 

Now the solution of (11.61) corresponding to Pj is 



and all eigenvalues need to be negative for the similarity solution to be stable. In 
that case, convergence to the fixed point is exponential, or algebraic in the original 
time variable t'. Soon the solution has effectively reached the fixed point, and there 
is very little change in the self-similar behaviour. If one or several of the eigenvalues 
around the fixed point vanish, the approach to the fixed point is slow, and the dynamics 
is effectively described by a dynamical system whose dimension corresponds to the 
number of vanishing eigenvalues. The same holds true if the attractor has few dimensions 
(such as a limit cycle or a low- dimensional chaotic attractor). Thus although singular 
behaviour is in principle a problem to be solved in infinite dimensions, in practise it 
typically reduces to a dynamical problem of few dimensions. In this review we analyse 
singularities from the point of view of the slow dynamics contained in (11.41) . to obtain an 
overview and tentative classification of possible scaling behaviours. We also emphasise 
the physical significance of these different types of behaviours. 

The perspective described above suggests a close relationship to the description of 
scaling phenomena by means of the renormalisation group, developed in the context of 
critical phenomena [23|, [21] ; we will continue to point out similarities, but we are not 



oo 





where uj is the eigenvalue: 
CP, = v,P,. 




p = e^^^Pi 
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aware that a classification similar to ours has been achieved using the language of the 
renormalisation group. For a computational perspective on analysing (11.41) in terms of 
its slow dynamics, see j25j. Finally, another approach sometimes associated with the 
classification of singularities is catastrophe theory [26] . However, as far as we are aware 
catastrophe theory only yields useful results if the problem can be mapped onto a low- 
dimensional geometrical problem, which can in turn be rephrased in terms of normal 
forms of polynomials. This has been shown to be the case for wave problems such as 
shock formation and wave breaking [27], as well as singularities of the eikonal equation 
[28] and related problems |29j . 

In this paper we discuss the following cases: 

(I) Stable fixed points (section [2]) 

In this case the fixed point is approached exponentially in the logarithmic variable 
r, so the dynamics is described by the self-similar law (II. 2p . This pure power-law 
behaviour is also known as type-I self-similarity [30] . 

(II) Centre manifold (section [3]) 

Here one or more of the eigenvalues around the fixed point are zero. As a result, 
the approach to the fixed point is only algebraic, leading to logarithmic corrections 
to scaling. This is called type-II self-similarity [30] ; it characterises cases where the 
blow-up rate is different from what is expected on the basis of a solution of the 
type dm). 
(HI) Travelling waves (section H]) 

Solutions of ( 11. ip converge to h = t'"0(^ + cr), which is a travelling wave solution 
of (ll.4p with propagation velocity c. 

(IV) Limit cycles (section [5]) 

Solutions have the form h = t'^'tp [C,, r] with ip being a periodic function of period T 
in T. This is known as "discrete self-similarity" [151 [31], since at times r„ = tq + hT, 
n integer, the solution looks like a self-similar one. 

(V) Strange attractors (section [6]) 

The dynamics on scale r are described by a nonlinear (low-dimensional) dynamical 
system, such as the Lorenz equation. 

(VI) Multiple singularities (section [Tj) 

Blow-up may occur at several points (xo,to) (or indeed in any set of positive 
measure), in which case the description ( 11. 4p is not useful. We also describe cases 
where (11.20 still applies, and blow-up occurs at a single point, but the underlying 
dynamics is really one of two singularities which merge at the singular time. 

This paper's aim is to assemble the body of knowledge on singularities of equations 
of the type (ll.ip that is available in both the mathematical and the applied community, 
and to categorise it according to the types given above. In addition to rigorous results 
we pay particular attention to various phenomenological aspects of singularities which 
are often crucial for their appearance in an experiment or a numerical simulation. For 
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Equation 


Type 


Dynamics 


Section 


Free surface flow 


ht + V- (/i"VA/i) ± V{hPVh) = 


1,11 


stable ? 


12.1.11 


+ {h'u)x = 


I 






/ \ /7O \ / 1 ^ /7"1\ 

p{ut + uux) = {h^Ux)x/h^ - [h ^)x 




stable 


12.1.11 


h,= [hKxl{l + hlYl^]^, 

K = l/{h{l + hlYl^)-hxx/{l + hlf/^ 


1 


stable 




ht + {hu)x = 0, ut + uux = hxxx 


I 


stable 


12.4.21 


^J{x-^,Y+a{x,t) 2a 


11 


Vr = —V^ 


13.2.11 


{h^\ + (/i^n),. = 


111 


stable 


m 


Geometric evolution equations 


ht = hj[l + hi) - 1/h 


11 


Ur = —U^ 


13.1.11 


^l^t = ^l^ss-{n-l){l-rs)/^ 


11 


Uj. = —u^ 


13.1.11 


React ion- diffusion equations 


Ut - Au = f{u) 


11 


Ur = —U^ 


13.1.21 


ut-V ■ (|m|™Vm) = vP 


11 


unknown 


13.1.21 


+ V ■ {pVS - Vp) = 0, p = -AS 


11 


Ur = —U^ 


13.2.21 


Nonlinear dispersive equations 


Ut + UUx = 


1 


stable 




#t + AV^+ l-^lPip = 


1,11 


Ur = —U^/v 
Vr = —UV 


13J 


Ut + U^Ux + Mxxx = 


11 


unknown 


13.3.11 




1 


unknown 


13.3.11 


ut = 2fv, vt = -2fu, ft = P 


IV 


circle 


m 


Choptuik equations 


I, IV 


limit cycle 


m 




1,11 


unknown 





Fluid equations 


Mt + (m • V)m = -Vp + An, V • M = 


1, IV ? 


unknown 




+ V)m = -Vp, V ■ M = 


1, IV ? 


unknown 




+ UUx + VUy = -Px + Uyy, Ux + Vy = 


1 


stable 






Table 1. A summary of PDE's discussed in this paper. The first column gives the 
PDE in question, the second the type of dynamics near the fixed point according to 
the classification enumerated above. In the case of attracting fixed-point dynamics, it 
is classed as "stable" , otherwise the equation governing the slow dynamics is given. 
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example, what are the observable implications of the convergence onto the self-similar 
form (11.21) being slow? In most cases, we rely on known examples from the literature, but 
the problem is almost always reformulated to conform with the formulation advocated 
above. However, some examples are entirely new, which we will indicate as appropriate. 
For each of the above categories, we will present at least one example in greater detail, 
so the analysis can be followed explicitely. A concise overview of the equations presented 
in this review is given in Table [H 



2. Stable fixed points 

A sub-classification into self-similarity of the first and second kind has been expounded 
in [32l |33l EH ES]. Self-similar solutions are of the first kind if (11.21) only solves (11.11) 
for one set of exponents a, (3; their values are fixed by either dimensional analysis or 
symmetry, and are thus rational. Solutions are of the second kind if solutions (11.21) exist 
locally for a continuous set of exponents however, in general these solutions are 
inconsistent with the boundary or initial conditions. Imposing these conditions leads to 
a non-linear eigenvalue problem, whose solution yields irrational exponents in general. 



2.1. Self- similarity of the first kind 

Our example, exhibiting self-similarity of the first kind [35], is that of a solid surface 
evolving under the action of surface diffusion. Namely, atoms migrate along the surface 
driven by gradients of chemical potential, see FiglS The resulting equations in the 
axisymmetric case, where the free surface is described by the local neck radius h{x,t), 
are [37] : 

h 



1 
h 



where 



1 h 



(2.1) 



is the mean curvature. In (I2.1I) . (I2.2I) . all lengths have been made dimensionless using an 
outer length scale R (such as the initial neck radius), and the time scale i?^/Z}4, where 
D4 is a forth-order diffusion constant. 

Physically, it is important to point out that (12. ip describes the evolution of the 
free surface at elevated temperatures, above the so-called roughening transition. This 
implies that the solid surface is smooth and does not exhibit facets, coming from the 
underlying crystal structure. Above the roughening transition, a continuum description 
is still possible [38]. The study of these models has lead to a number of interesting 
similarity solutions describing singular behaviour of the surface, such as grooves [39] or 
mounds [40ll4T]. 
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Figure 2. SEM images illustrating the pinch-ofF of a row of rectangular troughs in 
silicon (top) |36j . The bottom picture shows the same sample after 10 minutes of 
annealing at 1100°C. The troughs have pinched off to form a row of almost spherical 
voids. The dynamics is driven by surface diffusion. 



At a time t' ^ 1 away from breakup, dimensional analysis implies that i = t'^/^ is 
a local length scale. This suggests the similarity form 

h{x,t)=t'^'^H{x' /t'^'^), (2.3) 

and thus the exponents of (11. 2p are fixed by dimensional analysis, which is typical 
for self-similarity of the first kind. Of course, the result (12. 3p also follows when directly 
searching for a solution of (12.11) in the form of (II. 2p . In other cases, a unique set of local 
scaling exponents is determined by symmetry [12]. The similarity form of the PDE 
becomes 



- -iH 

4^ 



1 

n 



H 



X 



[2A) 



where k is the mean curvature of H. 

Solutions of (12. 4p have been studied extensively in To ensure matching to a 
time-independent outer solution, the leading order time dependence must drop out from 
(12. 3p . implying that 

e-±oo; (2.5) 
the general form of this matching condition for self-similar solutions of the form (II. 2p is 

e^±oo. (2.6) 
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H(^) 
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Figure 3. The approach to the self-similar profile for equation (|2.ip . The dashed line 
is the stable similarity solution H{^) as found from (|2.4p . The full lines are rescaled 
profiles found from the original dynamics p.ip at hm = 10~^,10~^, and hm = 10~^, 
respectively. As the singularity is approached, they converge rapidly onto the similarity 
solution (j^ . 



i g^(0) Q 

0.701595 1.03714 

1 0.636461 0.29866 

2 0.456842 0.18384 

3 0.404477 0.13489 

4 0.355884 0.10730 

5 0.326889 0.08942 



Table 2. A series of similarity solutions of (|2.4p as given in [43]. The higher-order 
solutions become successively thinner and flatter. 



All solutions of the similarity equation (12.11) . and which obey the growth condition (12.51) 
are symmetric, and form a discretely infinite set |13], similar to a number of other 
problems discussed below. The series of similarity solutions is conveniently ordered by 
descending values of the minimum, see table El Only the lowest order solution i?o(0 is 
stable, and is shown in Fig. [3l we return to the issue of stability in section 12.51 below. 
The fact that permissible similarity solutions form a discrete set implies a great deal 
of "universality" in the way pinching can occur. It means that the local solution is 
independent of the outer solution, and rather that the former imposes constraints on 
the latter; in particular, the prefactor c in (12.51) must be determined as part of the 
solution (see Table [2]). 
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2.1.1. Thin films and thin jets A further class of solutions displaying self-similarity of 
the first kind is the generalised long-wave thin-film equation 

ht + V ■ {h^'VAh - Bh'^Vh) = , n > 0. (2.7) 

The mass flux in this equation has two contributions: the first is due to surface tension, 
and the second is due to an external potential. When n = m = 3, then z = /i(x, t) 
represents the height of a film or a drop of viscous fluid over a flat surface, located at 
z = 0; the external potential is gravity. If B is negative, (12.71) describes a film that is 
hanging from a ceiling. Regardless of the sign of B, there is no singularity in this case 
|44j . The case n = 1 and 5 = corresponds to flow between two solid plates, to which 
we return in section 17.11 below. 

Solutions to (I2.7P are said to develop point singularities if h goes to zero in finite 
time. This happens if one incorporates van der Waals forces, which at leading order 
implies n = 3 and m = — 1 with B < 0. In [15], [16] (see also the review ^47j, where 
further full numerical simulations and mathematical theory are reported) the existence 
of radially symmetric self-similar touchdown solutions of the form 

h{r,t)=t'"^H{0, ^ = r/t'l (2.8) 

is shown numerically in this case. Self-similar solutions that touch down along a line 
exist as well, but they are unstable. A proof of formation of singularities in this context 
has been provided by Chou and Kwong [i8] . 

A related set of equations are those for thin films and jets, but which are isolated 
instead of being in contact with a solid. Problems of this sort furnish many examples 
of type-I scaling, as reviewed from a physical perspective in [19]. If the motion is no 
longer dampened by the presence of a solid, inertia often has to be taken into account. 
This means that a separate equation for the velocity is needed, which is essentially the 
Navier-Stokes equation below, but often simplified by a reduction to a single dimension. 
Thus one has solutions of the form 

h{x, t) = f'^HiO, u{x, t) = t'f'-'UiO, (2.9) 

where ^ = x'/t'^. If a > (3 the profile is slender, and the dynamics is well described in 
a shallow-water theory. In this case the equations for an axisymmetric jet with surface 
tension become 



dth'^ + d^{uh^) = (2.10) 



and 



8 (8 vh'^) 

p{8tu + u8,,u) = -{^/p)8,{l/h) + 3^^M^^ZLi. (2.11) 

The system (I2.10p . (l2.1ip is interesting because it exhibits different scaling 
behaviours depending on the balance between the three different terms in (12. lip [32]. 
This is an illustration of the principle of dominant balance, which is of great practical 
importance in practise, where it is a priori not known which physical effect will be 
dominant. In the case of (12. lip , these are the forces of inertia on the left, surface 
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tension (first term on the right), and viscosity (second term on the right). Pinching 
is driven by surface tension, so it must always be part of the balance. Three different 
possible balances remain [i2] : 

(i) In the first case [50], all forces in (12.111) are balanced as the singularity is 
approached. The exponents a = l,/3=l/2in (12.91) follow directly from this condition. 
As shown in [51], there is a discretely infinite sequence of self-similar profiles H{C,),U{^) 
corresponding to this balance. Numerical evidence strongly suggests that only the first 
profile, corresponding to the thickest thread, is stable [6]. All the other profiles are 
unstable, and thus cannot be observed. We will revisit this general scenario again 
below, when we study the stability of fixed points more generally. 

(ii) The second possibility corresponds to a balance between surface tension and 
viscous forces, thus putting p = in ( 12. lip . Physically, this occurs if the fiuid is very 
viscous [52j|. In section [2.4. II below we will describe the pinching solution corresponding 
to this case in more detail, as an example of self-similarity of the second kind. The 
exponent a = 1 is fixed by the balance, but f3 is fixed only by an integrability condition. 
This once more results in an infinite sequence of solutions, ordered by the value of f3. 
Again, only one profile, which has the largest value of /3 = 0.17487 is stable. This time, 
this corresponds to the smallest value of the minimum radius Rq, or the thinnest thread, 
as opposed to thickest thread in the case of the inertial-surface tension- viscous balance. 

If one inserts this viscous solution into the original equation (12.111) . one finds that 
in the limit t' —>■ 0, the inertial term on the left grows faster than the two terms 
on the right. This means that regardless how large the viscosity, eventually all three 
terms become of the same order, and one observes a crossover to the inertial-surface 
tension-viscous similarity solution described above, which is characterised by another 
set of scaling exponents and similarity profiles. In particular, the surface tension- viscous 
solution is symmetric about the pinch point, whereas the solution containing inertia is 
highly asymmetric [53] . We remark that crossover between different similarity solutions 
may also occur by another mechanism, not directly related to the dominant balance 
between different terms in the equation (cf. section [7!T]) . 

Equations (12. 101) . (12. Ill) correspond to a viscous liquid, surrounded by a gas, which 
is not dynamically active. The case of an external viscous fluid is considered in detail in 
section m below. The case of no internal fluid is special, in that the dynamics decouples 
completely into one for independent slices [53]. As a result, there is no universal proflle 
associated with the breakup of a bubble in a viscous environment, but rather it is 
determined by the initial conditions. 

(iii) At very low viscosity (z/ ^ in ( 12. lip ), the relevant balance is one where 
inertia is balanced by surface tension, so one might want to set z/ = in (12.111) . as 
done originally in [55]. However, the resulting equations do not lead to a selection of 
the values of the scaling exponents instead, there is a continuum of solutions [56], 
parameterised by the value of a, each with a continuum of possible similarity profiles. 
In fact, for vanishing viscosity ( l2.10l) .( l2rTT|) does not go toward a pinching solution, but 
the slope of the interface steepens, and one finds a shock solution ^57j, similar to the 
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generic scenario described in section 12.41 below. 

It was however shown numerically in |5H1 [59], and investigated in more detail 
in [60], that pinch-off of an inviscid fluid is well described by a solution of the full 
three-dimensional, axisymmetric potential flow equations. This is thus an example of 
a similarity solution of higher order in the independent variable, but both coordinate 
directions scale in the same way. The scaling exponents in (12.91) are a = jS = 2/3 in this 
case, which violates the assumption a > jS for the validity of the shallow water equations 
(12. lOp . ( 12.111) . In addition, we note that the similarity profile can no longer even be 
written as a graph as assumed in (12.91) . but turn over, as first observed experimentally 
in [61]. It is not known whether there also exists a sequence of similarity solutions, as 
in the case of the other balances. The case of no internal fluid is again very special, and 
leads to type-II scaling. It is considered in section 13.2.11 below. 

Finally, variations of (12. 101) . ( 12. Ill) have been investigated in [62]. Breakup was 
considered in arbitrary dimensions d (yet retaining axisymmetry) and with the pressure 
term 1/h replaced by an arbitrary power law l/h^. After introducing a new variable 
there remains a single parameter r = {d — l)/p, which can formally be varied 
continuously. For all values of r, discrete sequences of type-I solutions are obtained. 
For r > 1/2, profiles are asymmetric, while below that value they are symmetric. At 
the critical value, both types of solutions coexist. Another interesting feature of the 
limit r = 1/2 is that the viscous term becomes subdominant at leading order. However, 
similar to the case d = 3,p = 1 mentioned above, no selection takes place in the absence 
of the viscous term. Nevertheless, the solutions selected by the presence of the viscous 
term are very close to an appropriately chosen member of the family of inviscid solutions. 

2.2. Singularities in Euler and Navier-Stokes equations 

One of the most important open problems, both in physics and mathematics, is 
the existence of singularities in the equations of fluid mechanics: Euler and Navier- 
Stokes equations in three space dimensions. The Navier-Stokes equations represent the 
evolution of a viscous incompressible fluid and are of the form 

Ut + u-Vu = -Vp + Re-^Au, V • u = 0, (2.12) 

where u represents the velocity field, p the pressure in the fluid and Re is a dimensionless 
parameter called Reynolds number. Formally, by making Re— >^ oo, the term involving 
Au vanishes and we arrive at the Euler system, that models the evolution of the velocity 
and pressure fields of an inviscid incompressible fluid: 

Ut + u-Vu = -Vp, V-u = 0. (2.13) 

We exclude from our discussion certain "exact" blow-up solutions of the Euler equations 
[63], which have the defect that the velocity goes to infinity uniformly in space; in 
other words, they lack the crucial mechanism of focusing. Formally, they are of course 
similarity solutions of (12.131) . but with spatial exponent a = 0. 



Singularities 



13 



As we mentioned above, the existence of singular solutions is unknown. 
Nevertheless, some scenarios have been excluded. For the Navier-Stokes equations, 
there exists no nontrivial self- similar solution of the first kind 

u(x, t) = t'-^/^u (^) ^ ^ = x7t'^/2 (2.14) 

in L^(R^). This was proved by Necas, Ruzicka and Sverak [6l]. However, this does not 
exclude the formation of a singularity in a localised region: the matching condition (12.61) 
for this case implies |U| oc \^\~^ as |^| —>■ oo, which is not in L^. Therefore, the theorem 
[61] does not apply. 

A possible self-similar solution consisting of two skewed vortex-pairs has been 
proposed by Moffatt in [7] in the spirit of the scenario suggested by the numerical 
simulations of Pelz [65], of the implosion of six vortex pairs in a configuration with 
cubic symmetry. More recent numerical experiment by Hou and Li [66] seem to indicate 
that, although the velocity field may grow to very large values, singularities in the above 
mentioned scenarios saturate eventually and the solutions remain smooth. It has been 
argued in [67] that no self-similar solutions for Euler system should exist and that the 
"limit-cycle" scenario described in section [5] could apply. 

Under certain circumstances, such as special symmetry conditions or appropriate 
asymptotic limits, the Navier-Stokes and Euler systems may simplify and give rise to 
models for which the question of existence of singular solutions is somewhat simpler to 
analyse. This is the case for the Prandtl boundary-layer equations for the 2-D evolution 
of the velocity field (u, v) in y > 0: 

Ut + UUx + VUy = + Uyy, Ux+Vy = (2.15) 

with boundary conditions u = v = 0; p is a given pressure field and the behaviour of the 
velocity field at infinity is prescribed. Equation (12.151) describes the asymptotic limit 
of the Navier-Stokes equation near a solid body in the limit of large Reynolds numbers 
Re. The variable x measures the arclength along the body, and Re^^'^y is the distance 
from the body. Historically, a lot of attention was focused on the stationary version of 
(I2.15p . considering it as an evolution equation in x. At some position Xs along the body, 
the so-called Goldstein singularity t> oc (x^ — a;)~^/^ is encountered [68], which signals 
separation of the flow from the body. However, in reality the outer flow changes as a 
result of the appearance of a stagnation point, and one has to consider the interaction 
between the boundary layer and the outer flow [69] . 

It is thus conceptually simpler to consider the case of unsteady boundary layer 
separation, which is described by the first singularity of (12.151) at time to- The formation 
of singularities of f l2.15p in finite time was proved by E and Engquist [70]. It was first 
found numerically by van Dommelen and Shen \7T\, and its analytical structure was 
investigated in [72], using Lagrangian variables, which follow fluid particles as they 
separate from the surface (see also ^J). In the original Eulerian variables, the self- 
similar structure is [TH [75] 

, ,,1/2 ,1/2. s ^ x' - upt' y4>]^^ 
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where Mq, 4>o, and A are constants which depend on the problem, while U is universal 
and can be given in terms of elhptic integrals. Note that the exponents for u and x are 
the generic exponents for a developing shock (see section [23] below), while the similarity 
exponent in the ^/-direction is different from the scahng for two-dimensional breaking 
waves [27j. We stress that the appearance of a singularity in (12.151) does not mean that 
the full 2D Navier-Stokes equation has developed a singularity. Instead, lower order 
terms in the asymptotic expansion that lead to (12.151) become important close to the 
singularity. 

In relation with singularities in fluid mechanics, we can mention briefly a few 
important problems involving models or suitable approximations to the original Euler 
and Navier-Stokes systems. One concerns weak solutions to the Euler system for which 
the vorticity (c<j = V x u) is concentrated in curves or surfaces. This is the case of 
the so called vortex filaments and sheets in which the vorticity remains concentrated 
for all times, in absence of viscosity. A useful way to represent the vortex sheet, when 
it evolves in 2D, is by assuming the location of its points {x{a,t),y{a,t)) as complex 
numbers z{a,t) = x{a,t) + iy{a,t). Then, the evolution of z{a,t) is given by the 
so-called Birkhoff-Rott equation [76] : 

z*{a,t) = —PV r 7(y,^),^) ^ ^)^^/ (2.17) 
27ri z{a,t) - z{a',t) ' ^ ' ^ ' 

where z* stands for the complex conjugate of z. The principal value is denoted by PV, 
and 7 is the vortex strength and is such that dT = '-f{z{a,t),t)za{(y,t)da is constant 
along particle paths of the flow. The question then is whether or not these geometrical 
objects will remain smooth at all times or develop singularities in finite time. In the 
case of vortex sheets, singularities are known to develop in the form of a divergence 
of the curvature at some point. These are called Moore's singularities after their 
observation and description by D. W. Moore [77J. A mathematical proof of existence 
of these singularities is provided by Caflisch and Orellana in [78]. These singularities 
exhibit self-similarity of the first kind as shown, for instance, in [79|: if one defines the 
inclination angle 0{s,t) in terms of the arclength parameter s as such that Zg = e*^, 
then the curvature is given hj k, = Og and may blow-up in the self-similar form (up to 
multiplicative constants): 

K{s,t') = j^g{r]), v = s'/t', 0<6<1, (2.18) 



where 



oiv) = J- sin(5arctanr7) . (2-19) 

(1 + 772): 



2 



Interestingly, numerical simulations and Moore's original observations suggest that, 
although singular solutions with any 6 are possible, that the solution with 5 = | is 

preferred. Thus the generically observed geometry near the singularity is of the form 

3 

y = including the case of 3D simulations. This poses an interesting "selection 

problem" for the | power which has not received a definitive answer so far. 
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Another type of solution of fl2.17l) has the from of a double-branched spiral vortex 
sheet [80]. The explicit form is 

z(3t) = h"^' (2 20) 

where the two cases correspond to the two branches of the spiral. The parameter (3 
is related to integration variable a of (12.171) by d/S = Zada. The exponents are of the 
form u = 1/2 + ib and q = 1/2 + ifxb, corresponding to a vortex of radius r = t'^/^ 
collapsing in finite time. However, in this case the vortex sheet strength is found to 
increase exponentially at infinity [80] . 

Vortex filaments result as the limit of a vortex tube when the thickness tends to zero. 
The fluid flow around a vortex filament is frequently approximated by a truncation of the 
Biot-Savart integral for the velocity in terms of the vorticity. This leads to a geometric 
evolution equation for the filament (see [HI] , chapter 7, and references therein) that can 
be transformed, via Hasimoto transformation, into the cubic Nonlinear-Schrodinger in 
ID. This fact allowed Gutierrez, Rivas and Vega to construct exact self-similar solutions 
for infinite vortex filaments [HS]. One can also consider the vorticity concentrated in a 
region separating two fluids of different density and in the presence of gravitational 
forces. This is the case of the surface water waves system for which the existence of 
singularities is open [83] . 

A different approach in the study of singularities for Euler and Navier-Stokes 
equations in three space dimensions relies on the development of models that share 
some of the essential mathematical difficulties of the original systems, but in a lower 
space dimension. This is the case of the surface quasi-geostrophic equation popularised 
by Constantin, Majda and Tabak [84] : 

e^ + Y -Ve = 0, (2.21a) 

v = V^^, e = -{-Ay/^ip, (2.215) 

to be solved in c? = 2. This system of equations describes the convection of an active 
scalar 6, representing the temperature, with a velocity field which is an integral operator 
of the scalar itself. Nevertheless, the mere existence of singular solutions to this equation 
in the form of blow-up for the gradient of 6 is still an open problem. One-dimensional 
analogues of this problem, representing the convection of a scalar with a velocity field, 
which is the Hilbert transform of the scalar itself do have singularities in the form of 
cusps, as proved in [85], [86]. The structure of such singularities has been described in 
[87] and they are, in fact, of the type described in the next section, that is of the second 
kind. 



2.3. Self- similarity of the second kind 

In the example of the previous subsection, the exponents can be determined by 
dimensional analysis, or from considerations of symmetry, and therefore assume rational 
values. In many other problems, however, the scaling behaviour depends on external 
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parameters, set for example by the initial conditions. In that case, the scaling exponent 
can assume any value. Often, this value is fixed by a compatability condition, resulting 
in an irrational answer. We will call this situation self-similarity of the second kind 
[32| [35] . Since it is relatively rare that results are tractable analytically, we mention 
two simple examples for which this is possible, although they do not come from time- 
dependent problems. 

The first example is that of viscous flow near a solid corner of opening angle 2a 
[88] . For analogues of this problem in elasticity, see [HHl [90] as well as the discussion 
in [35]. This flow is described by a Stokes' equation, whose solution near the corner is 
expected to be 

i^ = r^fx{0). (2.22) 

If one of the boundaries is moving, scaling is of the first kind, and A = 2 (the so-called 
Taylor scraper pTj). However, if the flow is driven by two-dimensional stirring at a 
distance from the corner, A is determined by the transcendental equation 

sin2(A-l)a = -(A-l)sin2a. (2.23) 

If 2a < 146°, (12.231) admits complex solutions, which correspond to an infinite sequence 
of progressively smaller corner eddies. Since A is complex. The strength of the eddies 
decreases as one comes closer to the corner. 

The second example consists in calculating the electric field between two non- 
conducting spheres, where an external electric field is applied in the direction of 
the symmetry plane [92]. In this case the electric potential between the spheres is 
proportional to {p/{Rh))^~^, where p is the radial distance from the symmetry axis, 
R the sphere radius, and h the distance between the spheres. Thus in accordance with 
the the general ideas of self-similarity of the second kind, the singular behaviour is not 
controlled by the local quantity p/h, but the "outer" parameter R comes into play as 
well. We now explain two analytically tractable dynamical examples of self-similarity of 
the second kind. 



2.4- Breaking waves in conservation laws 

We only consider the simplest model for the formation of a shock wave in gas dynamics, 
which is Burger's equation 

ut + uu^ = 0. (2.24) 

It is generally believed that any system of conservation laws that exhibits blow up will 
locally behave like fl2.24p [94j|- For example. Fig. [H shows the steepening of a density 
wave in a gas, leading to a jump of the density in the picture on the right. In the words 
of [93]: "We conclude that an infinite slope in the theoretical solution corresponds to a 
shock in real life" . As throughout this review, we only consider the dynamics up to the 
singularity. Which structure emerges after the singularity depends on the regularisation 
used, as the continuation to times after the singularity is not unique [951 EE]. If the 
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Figure 4. Fringe pattern sliowing tlie steepening of a wave in a gas, leading to the 
formation of a shock, which is traveUing from left to right [93]. The vertical position of 
a given fringe is proportional to the density at that point. In the last picture a jump 
of seven fringes occurs. 



regularisation is diffusive, a shock wave forms [97]; if it is a tliird derivative, one finds a 
KDV soliton. Finally, regularisation by higher-order nonlinearities has been considered 
in [27] as a model of wave breaking. 

It is well known [98j that (12.241) can be solved exactly using the method of 
characteristics. This method consists in noting that the velocity remains constant along 
the characteristic curve 

z = uo{x)t + x, (2.25) 

where uo{x) = u{x,0) is the initial condition. Thus 

u{z,t)=uo{x) (2.26) 

is an exact solution to ( 12.24^ . given implicitly. 

It is geometrically obvious that whenever uo{x) has a negative slope, characteristics 
will cross in finite time and produce a discontinuity of the solution. This happens when 
dz/dx = 0, which will occur for the first time at the singularity time 

to = min(- ^ \ . 1, (2.27) 



dxUo{x 

at a spatial position x = Xm- This means a singularity will first form at 

Xo = x^- (2.28) 

Since (12.241) is invariant under any shift in velocity, we can assume without loss of 
generality that Uo{xm) = 0, and thus that Xq = Xm- This means the velocity is zero 
at the singularity. We now analyse the formation of the singularity using the local 
coordinates x',t'. In [27], this was done by expanding the initial condition uq in x', and 
using (I2.26p . using ideas from catastrophe theory [25]. Here instead we use the similarity 
ideas developed in this paper. 
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The local behaviour of fl2.24p near to can be obtained using the scaling 

u{x,t) = t"'U {x'/t'"+^) , (2.29) 

which solves (12.241) . The similarity equation becomes 

- aU + {I + a)^U^ + UU^ = 0, (2.30) 

with implicit solution 

e = -f/-Cf/^+^/". (2.31) 

The special case a = has the solution U = which is inconsistent with the matching 
condition (12. 6p . and thus has to be discarded. 

We are thus left with a continuum of possible scaling exponents a > 0, as is typical 
for self-similarity of the second kind. A discretely infinite sequence of exponents a„ 
is however selected by the requirement that ( I2.3ip defines a smooth function for all 
Namely, one must have 1 + 1/a odd, or 

«i = ^-^, ^ = 0,1,2..., (2.32) 

and we denote the corresponding similarity profile by Ui. The constant C in (12.311) 
must be positive, but is otherwise arbitrary. It is set by the initial conditions, which is 
another hallmark of self-similarity of the second kind. However, C can be normalised 
to 1 by rescaling x and U. We will see in section 1231 that the solution with ao, 

u{x,t)=t'^/^U^{x' If'/^) , (2.33) 

is the only stable one, all higher-order solutions are unstable. 

It is interesting to look at some possible exceptions to the form of blow-up given 
above, suggested by [M] : 

ut + uu^ = u". (2.34) 

This equation is also solved easily using characteristics. For a <2 the blow-up is alway 
of the form ( 12. 33^ . for cr > 2 two different behaviours are possible. For small initial 
data Mo(a;), a singularity still forms like (I2.33p . but in addition u may also go to infinity. 
However, there is a boundary between the two behaviours [91], where the slope blows 
up at the same time that u goes to infinity. For this case, one expects all terms in (I2.34p 
to be of the same order, giving 

M(x,t) =t'i^f/(0), ^ = x7t'^, (2.35) 
with similarity equation 

+ ^^^f/^ = f/" - UUi,. (2.36) 
i — a a — I 

The solution to (I2.36P that has the right decay at infinity is 
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Figure 5. The similarity solution (|2.37p for cr = 4. 



where C > is an arbitrary constant. The + and - signs describe the solution to the 

CT-2 

right and left of C,* = —{a — l)<'-i/(cr — 2), respectively. The special case a = 4 is 
shown in Fig. O The similarity solution fl2.37p is not smooth at its maximum; rather, 
its first derivative behaves like f/^ oc — This can be understood from the 

exact solution; in order for blow-up to occur at the same time that a shock is formed, 
the initial profile must already have a maximum with the same regularity as fl2.37p . 
Thus, the situation leading to fl2.35p is a very special one, requiring very peculiar initial 
conditions. 

2.4-1 ■ Viscous pinch-off As explained in section B.l.ll the pinch-off of a very viscous 
fluid is described by fl2.10p . (12. lip , with p = 0, but only for finite range of scales. 
The equations can be simplified considerably by introducing Lagrangian variables, i.e. 
writing all profiles as a function of a particle label s. This means the particle is at 
position z{s,t) at time t, and Zt{s,t) is the velocity at time t. The jet profile can be 
obtained from Zs = l/h'^{s,t), and (12. lip becomes 

The typical velocity scale is 'y/rj, where 7 is the surface tension and rj is the viscosity; 
(I2.38P has been made dimensionless accordingly. The time-dependent constant of 
integration C{t) has to be determined self-consistently. Note that the self-similar form 
(II. 2p is a solution of (I2.38P for a = 1, and any value of /3; the exponent /3 will be 
determined by the consistency condition (12.450 below. 
Since a = 1, a scaling solution of ( ]2.38p has the form 

h-\s, t) = t'-^f (0 , with ^ = s'/t'"' (2.39) 

and 



C{t) = -Cot' 



(2.40) 
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Figure 6. A drop of viscous fluid falling from a pipette 1 mm in diameter [99]. Note 
the long neck. 



The relationship with the exponent (3 defined in (12. 9p is simply /3 = 7 — 2, as found 
from passing from Lagrangian to Eulerian variables. Inserting f l2.39l) .( l2.40p into (12.381) 
we obtain 



Co, 



(2.41) 



where Cq is a constant. Imposing symmetry and regularity of /, we expand /(^) in the 
form 



MO = Rf + e^' + o{e^') , z = 0, 1, 2, . . . (2.42) 

where we have normalised the coefficient of ^2«+2 ^^-^^ This is consistent, since any 
solution of (I2.38P is only determined up to a scale factor. Instead, the axial scale is fixed 
by the initial conditions. The parameter Rq is the rescaled minimum of the profile: 
hra = Rot'. Inserting (12.421) into (I2.4ip . at order ^2*+^ one obtains 

1 „ 127-1 



Ro 



Co 



(2.43) 



12(7-1)' 24(7-1)2 

where we have put 7 = (i + 1)7. 

Each choice of i corresponds to one member in an infinite sequence of similarity 
solutions. Equation (I2.4ip can easily be integrated in terms of In ^ and y = a//: 



dy 



— Inf+^ + C 



iil + QRo)y^-y^-6Roy) QRo7 ' QRol 
with C an arbitrary constant. Computing the integral above we obtain 

y~^ ((27 -l)y + If"- {l-yr-= (2.44) 

which is an implicit equation for the i-th similarity profile y = yi[0) = \/ fiiO- 

The value of the velocity Uoo at infinity must be a constant to be consistent with 
boundary conditions. It can be found by integrating Zts = {h~^)t = t'~^(2/ + 7^/^) from 
zero to infinity: 



ztsds 



t'^~ 



24 (7-1)2 



(2.45) 
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where we have used f l2.41l) . The above condition Uoo = 0, which ensures that Uoo does 
not diverge as t' — > 0, is the equation which determines the exponent 7. Taking the 
derivative of (12.441) we obtain 

+ = ^ (r^ ((27 -i)y + If"' a-yV 

= -y-^~' (2^7 -y + If -77=^ 

v(i -y) 

which can be used to transform the integral in (12.451) to the variable y: 

3U^ 1 fWf^21-'^\.A 3 



" (l2(?3T))i = ITT X [[2—) ' - " 

/ i+1+7 1 2a-2~+3 1 2i+l \ 

(y-^+r- ((27 -l)y+ 1)^5^+^ (1 - — j dy = 0. (2.46) 
The function Ki{'y) may be written explicitly as 

, r(4-7)r(^) ^i(2^ + 2)7-r 



r(4-7+2lT2) V2 (^ + 1)7-1 

2^ + 3 1 ... , .r(3-7)r(2^) 



F{ 7 4-7; 4-7 + ;l-(2z + 2)7 -7 , , 

+ ^ 7,3-7;3-7 + 7T^;l-(2^ + 2)7l , (2.47) 



^2i + 2 " ' ' 2i + 2 
where F[a,b] c, z) is the hypergeometric function [100] . Roots of 7j are given in Table 

El 

To summarise, each exponent 7^ corresponds to a new member fi{C) of an infinite 
hierarchy of similarity profiles, to be found from (12.441) . If one converts the Lagrangian 
variables back to the original spatial variables, one obtains 

h{x,t) =t'(f)^st i^'/t'"'-^) ■ (2-48) 
Thus for t' ^ the typical radial scale t' of the generic i = solution rapidly becomes 
smaller than the axial scale t"^-^'^^ (cf. Table E])- This explains the long necks seen in 
Fig. El 

2.4-2. More examples Other recent examples for scaling of the second kind have been 
observed for the breakup of a two-dimensional sheet with surface tension. In a shallow- 
water approximation, which is justified for a description of breakup, the equations read 

m 

ht + {hu)x = 0, ut + uu^ = h^^^ (2.49) 

after appropriate rescaling. Local similarity solutions can be found in the form 

h{x, t) = t'^^-^H{r]), u{x, t) = t'f-^Uir]), (2.50) 

where rj = x' /t'^. The exponent (3 is not determined by dimensional analysis. Instead, 
it must be found from a solvability condition on the nonlinear system of equations for 
the similarity functions H, U. 
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i 


7i 


7") 





2.1748 


0.0709 


1 


2.0454 


0.0797 


2 


2.0194 


0.0817 


3 


2.0105 


0.0825 


4 


2.0065 


0.0828 


5 


2.0044 


0.0832 



Table 3. A list of exponents, found from Ki{'y) — using MAPLE, with Ki given by 
(|2.47p . The number 2i + 2 gives the smallest non- vanishing power in a series expansion 
of the corresponding similarity solution around the origin. Only the solution with i = 
is stable. The rescaled minimum radius is found from (j2.43[) . 



The result of the numerical calculation is |101] /3 = 0.6869 ± 0.0003, which is 
curiously close to f3 = 2/3, which is the value that had been conjectured earlier [102j . 
but contains a small correction. The value (3 = 2/3 comes out if both length scales in 
the longitudinal and transversal directions are assumed to be the same, implying that 
4/? — 2 = /?. This is a natural expectation for problems governed by Laplace's equation, 
such as inviscid, irrotational flow [59] , and indeed is observed for three-dimensional drop 
breakup [58l [59] . However, in present case, even if the full two-dimensional irrotational 
flow equations are used, /3 ^ 2/3. 

Other physical problems which frequently involve anomalous scaling exponents are 
strong explosions on one hand, and collapse of particles or gases into a singular state on 
the other. These types of problems have been reviewed in great detail in a number of 
textbooks and articles [32l [3H [33], [35] , but continue to attract a great deal of attention. 
As with many other singular problems, the type of scaling depends on the details of 
the underlying physics, and scaling of both the first and second kind is observed. For 
example, the radius of a shock wave resulting from a strong explosion can be calculated 
from dimensional analysis to be cx t^/^ |103j . However, in the seemingly analogous 
case of a strong implosion, an anomalous exponent is observed, which moreover depends 
on the parameters of the problem [104[ [98] . Cases were collapse and shock formation 
coincide were given by [105] (similar to section 12.4! above). In a somewhat different 
context, anomalous scaling is observed in model calculations for the collapse of self- 
gravitating particles [106] and Bose- Einstein condensates |107j . It is important to 
remember that these examples come from kinetic equations describing the stochastic 
collision of waves or particles, and hence involving nonlocal collision operators. However, 
the kinetic equations appear to be closely related to certain PDE problems pilSj , which 
are analogous to other evolution equations studied in this article. 
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2.5. Stability of fixed points 

Self-similar solutions correspond to fixed points of the dynamical system (11.4p . whose 
stability we now investigate by linearising around the fixed point. We explain the 
situation for the example of section 12.11 in more detail, for which the transformation 
reads 

hix,t)=t"/^H{^,T), (2.51) 
where r = — In(t'). The similarity form of (12. ip becomes 



H 



(2.52) 



which reduces to (12.41) if the left hand side is set to zero. To assure matching of (12.521) to 
the outer solution, we have to require that (I2.5ip is to leading order time- independent 
as ^ is large, which leads to the boundary condition 

Hr-{H - iH(:)/4: for |^| ^ oo. (2.53) 

This is the natural extension of (12. 5p to the time-dependent case. 

Next we linearise around any one of the similarity solutions H{S,) = Hi{^) listed 
in Table [2], as described in the Introduction. The stability is controlled by eigenvalues 
of the eigenvalue equation (ll.Sp . Inserting the eigensolution (ll.9p into (I2.53P one finds 
that Pj must grow at infinity like 

P.iO^e-'^^. (2.54) 

Similarly, the growth condition for the general case of a similarity solution of the form 
(O is 

P,(Ofxe^. (2.55) 

If the similarity solution H{$^) is to be stable, the real part of the eigenvalues of C 
must be negative. However, there are always two positive eigenvalues, which are related 
to the invariance of the equation of motion (12. ip under translations in space and time, 
as noted by [1091 IllOj . Namely, for any e, the translated similarity solution 

h^^\x,t)=t''/m{^) (2.56) 

is an equally good self-similar solution of (12. ip . and thus of (12.520 . In particular, we can 
expand (12.560 to lowest order in e, and find that 

/7(^)(e, r) = H{0 + ee^'H^iO + 0{e'), (2.57) 

where the linear term is a solution of (11.60 . 
Thus 

(e^^Fg)^ = e^^pHi: = e^^CH^. (2.58) 

But this means that = (3 = 1/4 is an eigenvalue of C with eigenfunction H^{^). 
Similarly, considering the transformation t — >■ t + e, one finds a second positive 
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eigenvalue z/j = 1, with eigenfunction ^H^. However, these two positive eigenvalues 
do not correspond to instability. Instead, the meaning of these eigenvalues is that upon 
perturbing the similarity solution, the singularity time as well as the position of the 
singularity will change. Thus if the coordinate system is not adjusted accordingly, it 
looks as if the solution would flow away from the fixed point. If, on the other hand, the 
solution is represented relative to the perturbed values of Xq and to, the eigenvalues 
and Ut will not appear. 

The eigenvalue problem (11. 8p was studied numerically in [13]. It was found that each 
similarity solution Hi has exactly 2i positive real eigenvalues, disregarding Vx-,i^t- The 
result is that the linearisation around the "ground state" solution Hq only has negative 
eigenvalues while all the other solutions have at least one other positive eigenvalue. This 
means that Hq is the only similarity solution that can be observed, all other solutions 
are unstable. Close to the fixed point, the approach to Hq will be dominated by the 
largest negative eigenvalue Vi. 

h{x,t)=t"/'[m) + ^t'~'''PiiO]- (2.59) 
For large arguments, the point C,cr where the correction becomes comparable to the 
similarity solution is ^ ~ et'~'^^^^~'^'^^ , and thus ~ t'~^^'^. This means that the region 
of validity of H{C,) expands in similarity variables, and is constant in real space. This 
rapid convergence is reflected by the numerical results reported in Fig. [31 More formally, 
one can say that for any e there is a 5 such that 

\h{x,t)-t'^/^H{0\<e (2.60) 

if < (5 uniformly as t' — > 0. 

We suspect that the situation described above is more general: the ground state is 
stable, while each following profile has a number of additional eigenvalues. In the case 
of the sequence of profiles Hi of (12.41) . two new positive eigenvalues appear for each new 
profile, corresponding to a symmetric and an antisymmetric eigenfunction. Below we 
give two more examples of the same scenario, for which we are able to give a simple 
geometrical interpretation for the appearance of two additional positive eigenvalues at 
each stage of the hierarchy of similarity solutions. The simplest case is that of shock 
wave formation (cf. section [23]) , for which everything can be worked out analytically. 

The dynamical system corresponding to the self-similar solution (12.291) is 

Ur - aU + {I + a) ^Ui: + UUi: = 0, (2.61) 
and so the eigenvalue equation for perturbations P around the base profile Ui becomes 

(«i - u)P - (1 + ai)^P^ - P(I7,)5 - P^Ui = 0, ^ = 0, 1, . . . (2.62) 
Here Ui is the ith similarity function defined by (I2.3ip for the exponents as given by 

The eigenvalue equation (12.621) is solved easily by transforming from the variable ^ 
to the variable U, using (12.311) : 

[a, + (2z + 3)I7f +') + ll = ^ Urn + (1 + a.pf''] , (2.63) 



P 
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with solution 

jj3+2i-2u(i+l) 

P = (2.64) 

l + {2^ + 3)uf^' 

The exponent 3 + 2i — 2i^(i + 1) must be an integer for (I2.64p to be regular at the origin, 
so the eigenvalues are 

2t + A-j_ .^12,... (2.65) 
^ 2z + 2 ^ ^ 

As usual, the eigensolutions are alternating between even and odd. However, we are 

interested in the first instance, given by (12.271) . at which a shock forms. This implies 

that the second derivative of the profile must vanish at the location of the shock, and 

the amplitude of the j = 3 perturbation must be exactly zero. 

Thus for i = the remaining eigenvalues are u = 3/2,1,0,-1/2,...; the first 
two are the eigenvalues = (3 = 1 + a and Vt = 1 found above. The vanishing 
eigenvalue occurs because there is a family of solutions parameterised by the coefficient 
C in ( 12.311) . All the other eigenvalues are negative, which shows that the similarity 
solution (12.331) is stable. In the same vein, for ai = 1/4 there are two more positive 
exponents: v = 5/4, 1, 1/2, 1/4, so the solution must be unstable. The same is of course 
true for all higher order solutions. Thus in conclusion the ground state solution Uq 
given by (I2.33P is the only observable form of shock formation. The same conclusion 
was reached in [27] by a stability analysis based on catastrophe theory. 

The sequence of profiles for viscous pinch-off, found in sect ion [2T3| suggests a simple 
mechanism for the fact that two new unstable directions appear with each new similarity 
profile of higher order. In fact, the argument is strikingly similar to that given for shock 
formation. Differentiating ( ]2.38p with respect to s one finds that a local minimum point 
Smin remains a minimum. Thus the local time evolution of the profile can be written as 

oo 

h{s, t) = h^ + Y^ Bj{t)s'K (2.66) 

For generic initial data -82(0) 7^ 0, so there is no reason why B2 should vanish at 
the singular time, which means that the self-similar solution /q will develop, which 
has a quadratic minimum. This situation is structurally stable, so one expects the 
eigenvalues of the linearisation to be negative. If however the coefficients Bj{0) are 
zero for j = 2, . . .2n — 1, they will remain zero for all times. Namely, if the first k 
s-derivatives of h vanish, one has 

^^'^* = -^' i = l---^' (2-67) 

so the first k derivatives will remain zero. Thus to find the similarity profile with i = I, 
one needs -62(0) = -83(0) = as an initial condition. This is a non-generic situation, 
and a slight perturbation will make B2 and B^ nonzero. In other words, there are two 
unstable directions, which take the solution away from fi{C,), as defined by (12.421) . In the 
general case, the linearisation around /j(0 "wiH have 2i positive eigenvalues (apart from 
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the trivial ones). Extensive numerical simulations of drop pinch-off in the inertial-surface 
tension-viscous regime (cf. section 12.1.11) suggests that the the hierarchy of similarity 
solutions again has similar properties in this well, although stabihty has not 

been studied theoretically. The ground-state profile is stable, while all the others are 
unstable [3^. Even when using a higher-order similarity solution as an initial condition, 
it is immediately destabilised, and converges onto the ground state solution fSl] . 

3. Centre manifold 

In section |2] we described the generic situation that the behaviour of a similarity solution 
is determined by the linearisation around it. In the case of a stable fixed point, 
convergence is exponentially fast, and the observed behaviour is essentially that of the 
fixed point. In this section, we describe a variety different cases where the the dynamics 
is slow. In all cases we are able to associate this slow dynamics with a fixed point in 
the appropriate variable(s), around which the eigenvalues vanish. Instead, higher-order 
non- linear terms have to be taken into account, and the slow approach to the fixed point 
is determined by a low- dimensional dynamical system. 
We consider essentially two different cases: 

(a) The dynamical system (11. 4p possesses a fixed point Ho{^), which has a vanishing 
eigenvalue, with corresponding eigenfunction The dynamics in the slow 
direction ip is described by a nonlinear equation for the amplitude a(r), which 
varies on a logarithmic time scale: 

h = t'" [HoiO + a{T)m] , C = x'/t'^. (3.1) 

(b) The dynamical system does not possess a fixed point, but has a solution of a slightly 
more general form: 

h = ho{T)H{0, ^ = x'/W{t), (3.2) 

where Hq and W are not necessarily power laws. To expand about a fixed point, 
we define the generalised exponents 

a = -dM/ho, (3 = -drW/W (3.3) 

which now depend on time. In the case of a type-I similarity solution, this 
reduces to the usual definition of the exponent. In the cases considered below, one 
derives a finite dimensional dynamical system for the exponents a,j3 (potentially 
including other, similarly defined scale factors). Once more, the exponents vary on 
a logarithmic time scale, which can be understood from the fact that the dynamical 
system possesses a fixed point with vanishing eigenvalues. 

Zero eigenvalues can also be associated to symmetries of the singularity, like 
rotational or translational symmetries, which lead to the existence of a continuum of 
similarity solutions. Another example, which concerns the dynamics inside the singular 
object itself, is wave steepening as described by fl2.3ip above. As seen from (12.651) . there 
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Figure 7. Nine images (of width 3.5 mm) sliowing liow a '^He crystal "flows" down 
from the upper part of a cryogenic cell into its lower part [112] . The recording takes a 
few minutes, the temperature is 0.32 K. 11 mK. The crystal first "drips" down, so that 
a crystalline "drop" forms at the bottom (a to c); then a second drop appears (d) and 
comes into contact with the first one (e); coalescence is observed (f) and subsequently 
breakup occurs (h). 



indeed is a vanishing eigenvalue associated with this continuum of solutions. Below 
we will not be concerned with this case, but only consider approach to the singularity 
starting from nonsingular solutions. 

3.1. Quadratic non-linearity: geometric evolution and reaction- diffusion equations 

The appearance of this type of nonlinearity is characteristic for various nonlinear 
parabolic equations and systems. The blow-up behaviour is characterised by the 
presence of logarithmic corrections in the similarity profiles. 

3.1.1. Geometric evolution equations: Mean curvature and Ricci flows Axisymmetric 
motion by mean curvature in three spatial dimensions is described by the equation 

"-(tt%-1)- 

where h{x,t) is the radius of the moving free surface. A very good physical realization 
of fl3.4p is the melting and freezing of a ^He crystal, driven by surface tension [111] , see 
Fig. [71 As before, the time scale t has been chosen such that the diffusion constant, 
which sets the rate of motion, is normalised to one. A possible boundary condition for 
the problem is that h{0, t) = h{L, t) = R, where R is some prescribed radius. For certain 
initial conditions h{x,0) = ho{x) the interface will become singular at some time to, at 
which h{xo,to) = and the curvature blows up. The moment of blow-up is shown in 
panel h of Fig. [71 for example. 
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Inserting the self-similar solution fll.2p into (13 ■4p . one finds a balance for a = (3 = 
1/2. The corresponding similarity equation is 

One solution of fl3.5p is the constant solution = y^. Another potential solution is 
one that grows linearly at infinity, to ensure matching onto a time-independent outer 
solution. However, it can be shown that no solution to (13. Sp . which also grows linearly 
at infinity, exists |113l I114j . Our analysis below follows the rigorous work in |30j . 
demonstrating type-II self-similarity. In addition, we now show how the description 
of the dynamical system can be carried out to arbitrary order. 

The relevant solution is thus the constant solution, but which of course does not 
match onto a time-independent outer solution. We thus write the solution as 

hix,t)=t"/^[V2 + g{^,T)'^, (3.6) 

with r = — \n{t') as usual. The equation for g is then 

^^ = ^-^ + rT^-257^' (3-^) 

which we solve by expanding into eigenf unctions of the linear part of the operator 

Cg = g-^g^/2 + g^^. (3.8) 
It is easily confirmed that 

CH2,{^/2) = z/,if2^(e/2), 2 = 0, 1, ... , (3.9) 

where Hn is the n-th Hermite polynomial p 



and z/j = 1 — i. Thus the first eigenvalue is z/q = 1, which corresponds to the positive 
eigenvalue Ut coming from the arbitrary choice of to- The other positive eigenvalue 
eigenvalue does not appear, since we have chosen to look at symmetric solutions, 
breaking translational invariance. However, the largest non-trivial eigenvalue ui is zero, 
and the linear part of (13.71) becomes 
Oct ■ 

-^ = (l-t)a,, 1 = 0,1,.... (3.11) 

Thus all perturbations with i > 1 decay, but to investigate the approach of the cylindrical 
solution, one must include nonlinear terms in the equation for ai. 
If we write 

oo 

^?(e,r) = ^a,(r)i/2.(e/2), (3.12) 

the equation for ai becomes 

^ = -23/2a2 + 0(aia,), (3.13) 
dr 
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whose solution is 

ai = l/(23/2r). (3.14) 

Thus instead of the expected exponential convergence onto the fixed point, the 
approach is only algebraic. Since all other eigenvalues are negative, the r-dependence 
of the is slaved by the dynamics of ai. Namely, as we will see below, aj = 0{t^^), so 
corrections to fl3.13p are of higher order. To summarise, the leading-order behaviour of 
(13. 4p is given by 

h{x, t) = 1/2 + aiiT)H2{0] , (3.15) 

as was proven by [30j . 

Now we compute the specific form of the higher-order corrections to (I3.15p . which 
have not been worked out explicitly before. If one linearises around (I3.14p . putting 
ai = af^ + ei, one finds 

= Ci + other terms. (3.16) 

dr T 

This means that the coefficient A of ei = A/t"^ remains undetermined, and a simple 
expansion of in powers of yields an indeterminate system. Instead, at quadratic 
order, a term of the form ei = A In r/r^ is needed. Fortunately, this is the only place in 
the system of nonlinear equations for where such an indeterminacy occurs. Thus all 
logarithmic dependencies can be traced, leading to the general ansatz 

(n) Si ^ ^ (In tY 

k=i+l 1=0 

where Si and Siki are coefficients to be determined. The index n is the order of the 
truncation. 

The coefficients can now be found recursively by considering terms of successively 
higher order in in the first equation: 

^ = - 24V2aia2 + 22a^ - 

272V2af - 19lV2a^ + 192a^a2 (3.18a) 

^ = -a2- V2/4:al + 6al-8V2aia2. (3.186) 
dr 

The next two orders will involve the next coefficient 03. From (I3.18ap and (]3.18&p . one 



first finds 6121 and 62, by considering 0(r ^) and 0(r ^), respectively. Then, at order 
0(r-("+i)) in the first equation, where n = 3, one finds all remaining coefficients 6iki 
in the expansion ( 13.17P up to A; = n. At each order in r~^, there is of course a series 
expansion in In r which determines all the coefficients. 
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We constructed a MAPLE program to compute all the coefficients up to arbitrarily 
high order (10th, say). Up to third order in the result is: 



, V2 17 In (r) \/2 73 \/2 
1 /4 — H ^4-^^ — -V + 



r 



16 r2 16 r3 

867 In (r) ^2 289 (In (r))^ ^2 



128 



02 



128 r3 



■1/32^ + 



5 V2 17 In (r) ^2 
16 ~ ~ 64 73 ' 



and thus h{x,t) becomes 



V2 + ai(r) (-2 + e) + a^ir) (l2 - 12^^ + ^4) 



(3.19a) 
(3.196) 

(3.20) 



from which one of course immediately finds the minimum. To second order, the result 
is 

1 3 + 171nr" 
27 8r2 



hr 



(2t 



nl/2 



1 



(3.21) 



First, the presence of logarithms implies that there is some dependence on initial 
conditions built into the description. The reason is that the argument inside the 
logarithm needs to be non-dimensionalised using some "external" time scale. More 
formally, any change in time scale t = t/to leads to an identical equation if also lengths 
are rescaled according to h = /i/y^. This leaves the pref actor in fl3.2ip invariant, 
but adds an arbitrary constant tq to r. This is illustrated by comparing to a numerical 
simulation of the mean curvature equation fl3.4p close to the point of breakup, see Fig. [Hi 
Namely, we subtract the analytical result fl3.2ip from the numerical solution hm/{2\/t') 
and multiply by r^. As seen in FigJHl the remainder is varying slowly over 12 decades 
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in t'. If the constant tq is adjusted, this small variation is seen to be consistent with the 
logarithmic dependence predicted by fl3.2ip . 

The second important point is that convergence in space is no longer uniform as 
implied by (12.601) for the case of type I self-similarity. Namely, to leading order the 
pinching solution is a cylinder. For this to be a good approximation, one has to require 
that the correction is small: ^"^/t <^ 1. Thus corrections become important beyond 
^cr ~ which, in view of the logarithmic growth of r, implies convergence in a constant 
region in similarity variables only. As shown in |lllj . the slow convergence toward 
the self-similar behaviour has important consequences for a comparison to experimental 
data. 

Mean curvature flow is also an example of a broader class of problems called 
generically "geometric evolution equations". These are evolution equations intended 
to gain topological insight by flowing geometrical objects (such as metric or curvature) 
towards easily recognisable objects such as constant or positive curvature manifolds. The 
most remarkable example is the so called Ricci flow, introduced in |115] . which is the 
essential tool in the recent proof of the geometrisation conjecture (including Poincare's 
conjecture consequence) by Grigori Perelman. 

Namely, Poincare's conjecture states that every simply connected closed 3- 
manifold is homeomorphic to the 3-sphere. Being homeomorphic means that both 
are topologically equivalent and can be transformed one into the other through 
continuous mappings. Such mappings can be obtained from the flow associated to 
an evolutionary PDE involving fundamental geometrical properties of the manifold. 
Thurston's geometrisation conjecture is a generalisation of Poincare's conjecture to 
general 3-manifolds and states that compact 3-manifolds can be decomposed into 
submanifolds that have basic geometric structures. 

Perelman sketched a proof of the full geometrisation conjecture in 2003 using 
Ricci flow with surgery [116] . Starting with an initial 3- manifold, one deforms it 
in time according to the solutions of the Ricci flow PDE (13.221) we consider below. 
Since the flow is continuous, the different manifolds obtained during the evolution 
will be homeomorphic to the initial one. The problem is in the fact that Ricci flow 
develops singularities in finite time, one of which we describe below. One would like 
to get over this difficulty by devising a mechanism of continuation of solutions beyond 
the singularity, making sure that such a mechanism controls the topological changes 
leading to a decomposition into submanifolds, whose structure is given by Thurston's 
geometrisation conjecture. Perelman obtained essential information on how singularities 
are hke, essentially three dimensional cylinders made out of spheres stretched out along a 
line, so that he could develop the correct continuation (also called "surgery" ) procedure 
and continue the flow up to a final stage consisting of the elementary geometrical objects 
in Thurston's conjecture. 

Ricci fiow is defined by the equation 



dgij 




(3.22) 



dt 
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for a Riemannian metric gij, where Rij is the Ricci curvature tensor. The Ricci 
tensor involves second derivatives of the curvature and terms that are quadratic in 
the curvature. Hence, there is the potential for singularity formation and singularities 
are, in fact, formed. As Perelman poses it, the most natural way to form a singularity 
in finite time is by pinching an almost round cylindrical neck. The structure of this kind 
of singularity has been studied in [117] . By writing the metric of a (n + l)-dimensional 
cylinder as 

g = ds^ + ip^gcan , (3.23) 

where gcan is the canonical metric of radius one in the n— sphere S"", tpi^s, t) is the radius 
of the hypersurface {s} x S*" at time t and s is the arclength parameter of the generatrix 
of the cylinder. 

The equation for if) then becomes 

(3.24) 

In |117j it is shown that for n > 1 the solution close to the singularity admits a 
representation that resembles the one obtained for mean curvature flow: 

22 (n — 1) 2t'v^ 

Namely, fl3.24p admits a constant solution u{C,,t) = 1, and the linearisation around it 
gives the same linear operator (13.81) as for mean curvature flow. Thus a pinching solution 
behaves as 

w(e, r) = 1 + a(r)/f2(e/2) + o(r-i), (3.26) 
where the equation for a is a,- = — 8a^, with solution a = l/(8r). 

3.1.2. Reaction- diffusion equations The semilinear parabolic equation 

Mj- Am- |m|P"^m = (3.27) 

is again closely related to the mean curvature flow problem (13. 4p . Namely, disregarding 
the higher order term in hx, (13.41) becomes 

ht = - I. (3.28) 

h 

Putting u = 1/h one finds 

ut = + - '^ul/u, (3.29) 

which is (13.271) in one space dimension and p = 3, once more neglecting higher-order 
non-linearities. As before, (13.271) has the exact blow-up solution 

M = (p- l)i^t'-^^. (3.30) 

If 1 < p < Pc = where d is the space dimension, then there are no other 
self-similar solutions to (13.271) [18]. and blow-up is of the form (13.301) (see [118] . [119] 
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and |12U] for a recent review). As in the case of mean curvature flow, corrections to 
fl3.30P are described by a slowly varying amplitude a: 

u = t'V(P-i)(p _ 1)1^ [1 _ aH2{^/2) + 0{1/t')] , ^ = x' /t"/\ (3.31) 

where a obeys the equation 

ar = -Apa^. (3.32) 

This result holds in 1 space dimension. In higher dimensions, one has to replace x by 
the distance to the blow-up set. 

This covers all range of exponents (larger than one, because otherwise there is 
no blow-up) in dimensions 1 and 2. The situation if p > is not so clear: if 
p > 1 + I then there are solutions that blow-up and "small" solutions that do not 
blow-up. Nevertheless, the construction of solutions as perturbations of constant self- 
similar solutions holds for any d and any p > 1. A simple generalisation of fl3.27p results 
from considering a nonlinear diffusion operator, 

Ut-V ■ {\urVu) = uP (3.33) 

and now the blow-up character depends on the two parameters m and p, see [121] . 

3.2. Cubic non-linearity: Cavity breakup and Chemotaxis 

More complex logarithmic corrections are possible if the linearisation around the fixed 
point leads to a zero eigenvalue and cubic nonlinearities. 

3.2.1. Cavity break-up As shown in jl22j . the equation for a slender cavity or bubble 



IS 



-L ^{x-iY + a{x,t) 2a' 

where a{x,t) = h'^{x,t) and h{x,t) is the radius of the bubble. Dots denote derivatives 
with respect to time t. The length L measures the total size of the bubble. If for 
the moment one disregards boundary conditions and looks for solutions to fl3.34p of 
cylindrical form, a{x,t) = ao(t), one can do the integral to find 

So In 3.35 

V «o / 2ao 

It is easy to show that an an asymptotic solution of fl3.35p is given by 
t' 

ao oc (3.36) 

corresponding to a power law with a small logarithmic correction. Indeed, initial theories 
of bubble pinch-off [1231 1124] treated the case of an approximately cylindrical cavity, 
which leads to the radial exponent a = 1/2, with logarithmic corrections. 

However both experiment |125] and simulation |122] show that the cylindrical 
solution is unstable; rather, the pinch region is rather localised, see Fig. [91 Therefore, 
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Figure 9. The pinch-ofF of an air bubble in water [125j . An initially smooth shape 
develops a localised pinch-point. 

0.60 
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t' 

Figure 10. A comparison of the exponent a between full numerical simulations of 
bubble pinch-off (solid line) and the leading order asymptotic theory (|3.43p (dashed 
line). 




it is not enough to treat the width of the cavity as a constant L; the width W is itself 
a time-dependent quantity. In |122j we show that to leading order the time evolution 
of the integral equation f l3.34p can be reduced to a set of ordinary differential equations 
for the minimum Oq of a{x,t), as well as its curvature Oq. 

Namely, the integral in (13 .34^ is dominated by a local contribution from the pinch 
region. To estimate this contribution, it is sufficient to expand the profile around the 
minimum at z = 0: a{x, t) = oq + (ao/2)2;^ + 0{z^). As in previous theories, the integral 
depends logarithmically on a, but the axial length scale is provided by the inverse 
curvature W = (2ao/aQ)^/^. Thus evaluating (13.341) at the minimum, one obtains jl22] 
to leading order 

aoln(4VrV«o) = al/{2ao), (3.37) 

which is a coupled equation for oq and W. Thus, a second equation is needed to close 
the system, which is obtained by evaluating the the second derivative of (I3.34p at the 
pinch point: 

The two coupled equations (I3.37p . (l3.38l) are most easily recast in terms of the time- 
dependent exponents 

2a = — S^ao/cio; 25 = —dra'^/aQ, (3.39) 
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where P = a — S, so are generalisations of the usual exponents in fll.2p . The 
exponent 6 characterises the time dependence of the aspect ratio W. Returning to the 
collapse (13.351) predicted for a constant solution, one finds that a = 1/2 and 6 = 0. 
In the spirit of the the previous subsection, this is the fixed point corresponding to 
the cylindrical solution. Now we expand the values of a and 6 around their expected 
asymptotic values 1/2 and 0: 

a = l/2 + M(r), 5 = v{t). (3.40) 

and put w{t) = l/ln(aQ). 

To leading order, the resulting equations are 

Ut- = u + w/4:, Vr = —v — w/4:, Wr = 2vw'^. (3-41) 

The linearisation around the fixed point thus has the eigenvalues and —1, in addition 
to the eigenvalue 1 coming from time translation. As before, the vanishing eigenvalue 
is the origin of the slow approach to the fixed point observed for the present problem. 
The derivatives Ur and Vr are of lower order in the first two equations of (I3.4ip . and 
thus to leading order u = v and v = —w/i. Using this, the last equation of (13.411) can 
be simplified to 

= -w^/2. (3.42) 

Equation ( 13.42^ is analogous to ( 13.13^ . but has a degeneracy of third order, rather 
than second order. Equation (13.421) yields, in an expansion for small S [122] . 

a = 1/2 + ^ + 0(r), 6 = ^ + 0{r~y'). (3.43) 

Thus the exponents converge toward their asymptotic values a = (3 = 1/2 only very 
slowly, as illustrated in Fig. [TOl This explains why typical experimental values are 
found in the range a ~ 0.54 — 0.58 |125] . and why there is a weak dependence on initial 
conditions |126j . 

3.2.2. Keller-Segel model for chemotaxis This model describes the aggregation of 
microorganisms driven by chemotactic stimuli. The problem has biological meaning 
in 2 space dimensions. If we describe the density of individuals by u{x,t) and the 
concentration of the chemotactic agent by v{x,t), then the Keller-Segel system reads 

ut = Am - xV ■ (mVw), (3.44a) 
Tvt = Av + {u-l), (3.446) 

where F and x ^I's positive constants. In [13l I127j it was shown that for radially 
symmetric solutions of (I3.44ap .( r3.44&l) singularities are such that to leading order u 
blows up in the form of a delta function. The profile close to the singularity is self- 
similar and of the form 

^i^^t) = ( , (3.45) 



i?2(t) \R(t) 
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where 



and 



R{t) = Ce~^"--"'-2''^^+4^(i + o(l)) (3.46) 



^(0 - ^(TTa- 

The result comes from a careful matched asymptotics analysis that, in our notation, 
amounts to introducing the time-dependent exponent 

7 = -drR/R, (3.48) 

which has the fixed point 7 = 1/2. Corrections are of the form 

7=2 + f + (3-49) 

where a is controlled by a third-order non-linearity, as in the bubble problem: 

ar = -a^{l-a + o{a)). (3.50) 

3.3. Beyond all orders: The nonlinear Schrddinger equation 
The cubic nonlinear Schrodinger equation 

iipt + A(p + \(pf if = (3.51) 

appears in the description of beam focusing in a nonlinear optical medium, for which 
the space dimension is = 2. Equation (13.511) belongs to the more general family of 
nonlinear Schrodinger equations of the form 

i^t + Aip + ^ = 0, (3.52) 

and in any dimension d. Of particular interest, from the point of view of singularities, 
is the critical case p = A/d. In this case, singularities with slowly converging similarity 
exponents appear due to the presence of zero eigenvalues. We will describe this situation 
below, based on the formal construction of Zakharov |128j . later proved rigorously by 
Galina Perelman [129j . At the moment, the explicit construction has only been given 
for (i = 1, that is, for the quintic Schrodinger equation. The same blow-up estimates 
have been shown to hold for any space dimension c? < 6 by Merle and Raphael [130j . 
[131] . without making use of Zakharov's |128j formal construction. Merle and Raphael 
also show that the stable solutions to be described below are in fact global attractors. 
In the critical case (13.521) becomes in d=l: 

i(pt + ipxx + \<-pf^P = 0. (3.53) 

This equation has explicit self-similar solutions (in the sense that rescaling x Ax, 
t A^t, X^Lf leaves the solutions unchanged except for the trivial phase factor 

g-2winA^ of the form 



ifix, t) = e^/^°-e-^^^v.o(0, e = x'/t"/'. (3.54) 
t i 



Singularities 



37 



The function (po{^) solves 

- V'o,?? + V^o - Iv^ol'^V^o = 0, (3.55) 
and is given explicitly by 

^„({) = 'ff°'' . (3.56) 

cosh 2 (2y7lo';) 



We seek solutions of (13.531) using a generalisation of (13.541) . which allow for a 
variation of the phase factors, and the amplitude to be different from a power law: 

^{x, t) = e''^(*)-*''W"'/^A5 (t)<^,(2), (3.57) 

where z = X{t)x and ipa satisfies 

- V9a,^^ + V^a - ^az'^ifa - |v?a|^ = 0. (3.58) 

When h (= ^/a) is constant, (I3.57P is a solution of (I3.53P if (/i, A,/3) satisfy 

Ht = (3.59a) 

X-^Xt = f3 (3.595) 
+ X^f3^ = X^h\ (3.59c) 

Notice that the equation for /x is uncoupled, so we only need to solve the equations 
for (A,/5) simultaneously and then integrate the equation for yU. It is interesting for the 
following that, in addition to the solutions for constant a, one can let a vary slowly in 
time. The resulting system for (A, (3, h) is 

A-^Aj = p (3.60a) 
+ X^f3^ = X^h^ (3.605) 

ht = -cAV*°/V^- (3.60c) 

Note the appearance of the factor e~^°^'^ in the last equation, which comes from a 
semi classical limit of a linear Schrodinger equation with appropriate potential (see |129j ). 
and 

So= f ^Jl-s^/Ads=l. (3.61) 
Jo ^ 

Sq is an It follows from the presence of this factor that the non-linearity is beyond all 
orders, smaller than any given power, in contrast to the examples given above. 
As in section I3.2.H we rewrite the equations in terms of similarity exponents, 

to obtain the system: 

ttr = —(1 + 2a; + 7)0; (3.63a) 
7^ = (1 + 2a + 7)a - (7 + a)(l + 2a + 25 - 7) (3.635) 

6r = (-1 - 2a + 26)5 - (3.63c) 

h 

hr=- Sh. (3.63d) 
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The advantage of this formulation is that the exponents have fixed points. There are 
two famihes of equihbrium points for fl3.63flP - fl3.63djl : 

(1) a = — |, 7 = 0,5 = 0, h arbitrary positive or zero. 

(2) a = —1, 7 = 1,5 = 0, h arbitrary positive or zero. 

We first investigate case (1) by writing 

a = ~ + ai, 7 = 7i, S = 6i, h = hi. (3.64) 

The final fixed point corresponding to the singularity is going to be ai = 71 = 5i = 
hi = 0. However, there are also equilibrium points for any /i > 0, in which case the 
linearisation reads: 

ai,T = «! + -71 (3.65a) 

7i,r = -71 + ^1 (3.655) 

5^^^ = 25j - 2ai5i - 6fj-. (3.65c) 

This system has the matrix 

A 



1 


1 


1 

2 





\ 







-1 


1 




\ 











/ 



whose eigenvalues are: 1,0, and —1. The vanishing eigenvalue corresponds to the line 
of equilibrium points for /i > 0, the positive eigenvalue to the direction of instability 
generated by a change in blow-up time. The eigenvector corresponding to the negative 
eigenvalue gives the direction of the stable manifold. 

At the point h = 0, there is an additional vanishing eigenvalue, and the equations 
become: 

ai,T' = ("1 + ^7i)^i (3.66a) 
7i,r' = i-li + Si)hi (3.666) 
6iy = {26l - 2ai6i)hi - d^So (3.66c) 
hiy = — 5ihl, (3.66(i) 

where dr' = dr/hi. The first two equations reduce to leading order to 71 = 6ihi and 
ai = —6ih1/2, while the last two equations reduce to the nonlinear system: 

Siy = —S^Sq, hiy = —5ihl, Tr' = hi. (3.67) 

In the original r-variable, the dynamical system is 

5i,T = —SfSo/hi hiy = —5ihi, (3.68) 

which controls the approach to the fixed point. The system (13.681) is two-dimensional, 
corresponding to the two vanishing eigenvalues. 
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Integrating the first equation of fl3.67l) one gets Si ~ 1/{Sot'), and tlius using tlie 
second equation hi ~ 6*0/ In r'. From tlie last equation one obtains to leading order 
t' ~ rlnr/So, so tliat 

hi ~ — ^ , 5i ~ — — . (3.69) 
In r r In r 

Tlius we can conclude tliat 

air) ~ \ — , 7(r) ~ — , 6(t) ~ — . (3.70) 

In this fashion, one can construct a singular solution such that 

t .\ -jrlnT-ax'2/4(lnT)3 / (Iut)^ 

<l^^,(<!^il)ix) (371) 

Note the remarkable smallness of this correction to the "natural" scaling exponent of 
t'^/^, which enters only as the logarithm of logarithmic time r. 

The fixed points (2) can be analysed in a similar fashion. The linearisation leads 

to 

ai^r = 2ai + 7i (3.72a) 

7i,r = 71 (3.726) 

Si^r = Si. (3.72 c) 

All eigenvalues are positive, so one cannot expect these equilibrium points to be stable. 

One may also consider the blow-up of vortex solutions to both critical and 
supercritical solutions to nonlinear Schrdinger equation in 2D. These are a subset of 
the general solutions to NLSE that present a phase singularity at a given point. The 
singularities appear in the form of collapse of rings at that point. Both the existence of 
such solutions and their stability have been considered recently in |132[ 1133] . 

3.3.1. Other nonlinear dispersive equations The nonlinear Schrodinger equation 
belongs to the broader class of nonlinear dispersive equations, for which many questions 
concerning existence and qualitative properties of singular solutions are still open. 
Nevertheless, there have been recent developments that we describe next. 
The Korteweg-de Vries (KdV) equation 

ut + {u^x + u^)x = (3.73) 

describes the propagation of waves with large wave-length in a dispersive medium. For 
example, this is the case of water waves in the shallow water approximation, where 
u represents the height of the wave. In the case of an arbitrary exponent of the 
nonlinearity, (13.731) becomes the generalised Korteweg de Vries equation: 

Ut + {uxx + u'')x = , p > 1. (3.74) 
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Based on numerical simulations, [134j conjectured the existence of singular solutions 
of (13.741) with type-I self-similarity if p > 5. In [135j . |136j it was shown that in the 
critical case p = 5 solutions may blow-up both in finite and in infinite time. Lower 
bounds on the blow-up rate were obtained, but they exclude blow-up in the self-similar 
manner proposed by |134j . 

The Camassa-Holm equation 

Ut — Uxxt + SUxU = 2UxUxx + UxxxU (3.75) 

also represents unidirectional propagation of surface waves on a shallow layer of water. 
It's main advantage with respect to KdV is the existence of singularities representing 
breaking waves [137] . The structure of these singularities in terms of similarity variables 
has not been addressed to our knowledge. 



4. Travelling wave 

The pinching of a liquid thread in the presence of an external fluid is described by 
the Stokes equation |138j . For simplicity, we consider the case that the viscosity rj of 
the fluid in the drop and that of the external fluid are the same. An experimental 
photograph of this situation is shown in Fig. [H To further simplify the problem, we 
make the assumption (the full problem is completely analogous) that the fluid thread 
is slender. Then the equations given in [5] simplify to 

ht = -Vxh/2-vhx, (4.1) 



where 



1 / h\y) 



Wh\y) + {x-yy 



K dy, (4.2) 



y 



and the mean curvature is given by (12. 2p . Here we have written the velocity in units of 
the capillary speed = 7/77. The limits of integration X- and x+ are for example the 
positions of the plates which hold a liquid bridge 1139] . 

Dimensionally, one would once more expect a local solution of the form 

h{x,t)=t'H(^], (4.3) 



J' ^ 

and H{C,) has to be a linear function at infinity to match to a time- independent outer 
solution. In similarity variables, (14. 2 p has the form 

We have chosen real-space variable close to the pinch-point, such that the 

similarity description is valid in [—Xb-,Xb\. But if H is linear, the integral in ( 14. 4p 
diverges like 61nt', where 

1 H-i- H , . 

"=-4 irm^TTw ■ f^-^' 
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Here if+ and H are the slopes of the similarity profile at ±00. But this means that a 
simple "fixed point" solution (14.31) is impossible. 

However by subtracting the singularity as t' 0, one can define a self-similar 
velocity profile according to 

where now 

V{0 = V^^''\0-br, (4.7) 
and an arbitrary constant has been absorbed into V^^^\ In terms of V^^^\ and putting 

h{x,t)=t'H{^,T), (4.8) 
the dynamical system for H becomes 

H^ = H~{^ + V'^^''^) - W/^°V2 + brH^. (4.9) 
This equation has a solution in the form of a travelling wave: 

Hi^,T) = HiC). V^''^\^,T) = ViC), where C = ^ - br. (4.10) 
The profiles H, V of the travelling wave obey the equation 

H-{C + V)Hc = HVj2. (4.11) 
The numerical solution of the integro-differential equation (14.111) gives 

hmin = aontVr,t', whcrc a^nt = 0.033. (4.12) 
The slope of the solution away from the pinch-point are given by 

H+ = 6.6 and = -0.074, (4.13) 

which means the solution is very asymmetric, as confirmed directly from Fig. [1] These 
results are reasonably close to the exact result, based on a full solution of the Stokes 
equation [5]; in particular, the normalised minimum radius is flout = 0.0335 for the full 
problem. 



5. Limit cycles 

An example for this kind of blow-up was introduced into the literature in [15] in the 
context of cosmology. There is considerable numerical evidence |140] that discrete self- 
similarity occurs at the mass threshold for the formation of a black hole. The same 
type of self-similarity has also been proposed for singularities of the Euler equation 
|14H [67j . the porous medium equation driven by buoyancy [141], and for a variety of 
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other phenomena |142j . A reformulation of the original cosmological problem leads to 
the following system: 



fx 



1)/ 



X 



X ' 
(f + x)U' - if - x)V' 

X 

f[il-a^)U + V]-xUt 

x{f + x) 
f[il-a^)U + V]+xVt 



+ 1 



/a' 



(5.1a) 
(5.1&) 
(5.1c) 
{5.1d) 
(5.1e) 



x{f-x) 

In [in], the self-similar description corresponding to the system fl5.1al) - fl5.1ep was solved 
using formal asymptotics and numerical shooting procedures. This leads to the solutions 
observed in [15]. We now propose another system, which shares some of the structure 
of fl5.1ap -( !5Tej) . but which we are able to solve analytically: 

Ut{x,t) = 2f{x,t)v{x,t), (5.2a) 
vtix, t) = -2f{x, t)u{x, t), (5.25) 
ftix,t) = f{x,t). (5.2c) 

The system ( 15.2ap -( !5l2cj) is driven by the simplest type of blow-up equation ( ]5.2cp . and 
can be solved using characteristics. However, in the spirit of this review, we transform 
to similarity variables according to: 

w = f/(e,r) 



f = t'-'Fi^,r) 



(5.3a) 
{5.3b) 
(5.3c) 
> 0. Near 



It is seen directly from fl5.2cp that / first blows up at a local maximum /„ 
a maximum, the horizontal scale is the square root of the vertical scale t', and thus we 
must have ^ = x' jt'^l'^ . With that, the similarity equations become 

= + FV (5.4a) 

Vr = - FU {5Ab) 

Fr = -F~ iF^/2 + (5.4c) 

The fixed point solution of the last equation is 

^ (5.5) 

(5.6) 
(5.7) 



where c > is a constant. The equations for f/, V are solved by the ansatz 
U = Uo sin (C(0 + r) , V = Uo cos (C(0 + r) , 

and for the function C(^) one finds 
eC"(e)/2 = F - 1, 
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with solution C(^) = — ln(l + c^^). Thus (a single component of) the singular solution 
is indeed of the general form 

U = ^m)+r), (5.8) 

where ip is periodic in r. This is a particularly simple version of discretely self-similar 
behaviour, i.e. when T is the period of ip, the same self-similar picture is obtained for 
r = To + nT. 



6. Strange attractors and exotic behaviour 

In connection to limit cycles and in the context of singularities in relativity, a few 
interesting situations have been found numerically quite recently. One of them is the 
existence of Hopf bifurcations where a self-similar solution (a stable fixed point) is 
transformed into a discrete self-similar solution (limit cycle) as a certain parameter 
varies (see |143] ). Other kinds of bifurcations, for example of the Shilnikov type, are 
found as well [144] . Before coming to simple explicit examples, we mention that possible 
complex dynamics in r has long been suggested for simplified versions of the inviscid 
Euler equations |145i 11461 1141j . For a critical discussion of this work, see |147[ [8T]. 

The problems considered in these papers were the 2D axisymmetric Euler equations 
with swirl, which produces a centripetal force. In the limit that the rotation is confined 
to a small annulus, the direction of acceleration is locally uniform, and the equation 
reduces to that of 2D Boussinesq convection, where the centripetal force is replaced by 
a "gravity" force. Another related model is 2D porous medium convection, for which 
the equation reads 

^ + (Te,-V0)-VT = O, (6.1) 

where v = Te^ — V0 plays the role of the velocity field and T is the temperature. 
The potential follows from the constraint of incompressibility, which gives A(j) = Ty. 
Simulations provide evidence of a self-similar dynamics of the form |141] 

T = t'W(x7t'^+^r), (6.2) 

where r/ is approximately 0.1 and M is a function that is slowly varying with r. 

Depending on the model, both periodic behaviour as well as more complicated, 
chaotic motion has been observed in numerical simulations. Oscillations of temperature 
in r are motivated by the observation that a sharp, curved interface (i.e. the transition 
region between a rising "bubble" of hot fluid and its surroundings) becomes unstable and 
rolls up. However, owing to incompressibility, the sheet is also stretched, which stabilises 
the interface, leading to an eventual decrease in gradients. Locality suggests that this 
process could repeat itself periodically on smaller and smaller scales |141] . However, 
simulations of the Euler equation have also shown examples of a more complicated 
dependence on r, which might be chaotic behaviour |145] . We also mention that 
corresponding chaotic behaviour has been proposed for the description of spin glasses 
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in the theory of critical phenomena |148j . We now give some exphcit examples of chaos 
in the description of a singularity. 

In section 13.1.11 we treated a system of an infinite number of ordinary differential 
equations for the coefficients of the expansion of an arbitrary perturbation to an explicit 
solution. Such high-dimensional systems in principle allow for a rich variety of dynamical 
behaviours, including those found in classical finite dimensional dynamical systems, such 
as chaos. Consider for instance an equation for the perturbation g (the analogue of (13.71) ) 
of the form 

g, = Lg + F{g,g), (6.3) 

where Lg is a linear operator. Assuming an appropriate non-linear structure for the 
function F, an arbitrary nonlinear (chaotic) dynamics can be added. 

To give an explicit example of a system of PDE's exhibiting chaotic dynamics, 
consider the structure of the example given in section [5l It can be generalised to produce 
any low-dimensional dynamics near the singularity, as follows by considering the system 

M«(x,t) =2/F,({m«}), z = l,...,n, (6.4a) 
ftix,t) = fix,t). (6.46) 

Using the ansatz analogous to (15.61) : 

u» = [/W(C(0 + r,0, (6.5) 
and choosing C{^) = — ln(l + c^^), one obtains the system 

f/» = F, {[/«}. (6.6) 
To be specific, we consider n = 3 and 

Fi = a(n(2)-««), F2 = puW-M^^)-^^^^^^), F, = u'-'^^'^ -/3u^'\ (6.7) 

so that (16.61) becomes the Lorenz system [149] . As before, for t' 0, the variable 
r goes to infinity, and near the singularity one is exploring the long-time behaviour 
of the dynamical system (16.51) . In the case of (16. 7p . and for sufficiently large p, the 
resulting dynamics will be chaotic. Specifically, taking a = 10, p = 28, and j3 = 8/3, as 
done by Lorenz [150] . the maximal Lyapunov exponent is 0.906. The initial conditions 
with which ( 16. 5 p is to be solved depend on ^. Thus the chaotic dynamics will follow a 
completely different trajectory for each space point. As a result, it will be very difficult 
to detect self-similar behaviour of this type as such, even if data arbitrarily close to the 
singularity time is taken. If for example a rescaled spatial picture is observed at constant 
intervals of logarithmic time r, the spatial structure of the singularity will appear to 
be very different. However, as pointed out in [145] , chaotic motion is characterised by 
unstable periodic orbits, for which one could search numerically. 
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7. Multiple singularities 

The singularities described so far occur at a single point xq at a given time to- This need 
not be the case, but blow-up may instead occur on sets of varying complexity, including 
sets of finite measure. We begin with a case where singularity formation involves two 
different points in space. 

7.1. Hele-Shaw equation 

A particularly rich singularity structure is found for a special case of (12. 7p in one space 
dimension with n = 1. Dropping the second term on the right, which will typically be 
small, one arrives at 

ht + {hh^^^)^ = 0. (7.1) 

This is a simplified model for a neck of liquid of width h confined between two parallel 
plates, a so-called Hele-Shaw cell, which is a simplified model for the free surface in a 
so-called Hele-Shaw cell |151j . Breakup of a fluid neck inside the cell corresponds to h 
going to zero in finite time. 

Singular solutions displaying type-I self-similarity would be of the form 

/i(x,t) =t'°iy(x7t'("+^)/^), (7.2) 

but are never observed. Instead, several types of pinch solutions different from (17.21) 
have been found for (17. ip using a combination of numerics and asymptotic arguments 
[152[ I102[ 1153] . On one hand, singularities exhibit type-ll self-similarity. On the other 
hand, the simple structure (17. 2p is broken by the fact that the location of the pinch point 
is moving in space. The root for this behaviour lies in the fact that two singularities 
are interacting over a distance much larger than their own spatial extend. Below we 
report on three different kinds of singularities whose existence has been confirmed by 
numerical simulation of (17. ip . 

The first kind of singularity was called the imploding singularity in |153j . since it 
consists of two self-similar solutions which form mirror images, and which collide at the 
singular time. Locally, the solution can be written 

h{x, t) = t'^H{{x' + at')/t'^), (7.3) 

where —a is the constant speed of the singular point. Note that the scaling exponents 
do not agree with (17. 2p . The reason is that the singularity is moving, so h is the solution 
of 

hh^xx = Jit') = t'\ (7.4) 

where J is determined by matching to an outer region. The similarity profile if is a 
solution of the equation HH'" = 1, with boundary conditions 

Hij]) (X7f/2, 7] ^ -oo; H{7]) oc ^/8/3{A-r])^/^, T] ^ oo. (7.5) 

One might wonder whether this behaviour is generic, in the sense that it might 
depend on the initial conditions being exactly symmetric around the eventual point of 
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log (h) 




Figure 12. Same as FiglTl] but both parameters w — 0.07 and S ~ 0.01. 



blow up. The simulation of (17. ip shown in Fig. [TT] shows that this is not the case. The 
initial condition is 



h{x,0) = l-{l-w) 



3 6 1 

- cos Tcx COS 2iTx H COS 37rx(l + 6 sin 2ttx) 

2 10 10 ^ ' 



,(7. 



which for S = reduces to the symmetric initial condition considered by [153] . The 
type of singularity that is observed (or no singularity at all) depends on the parameter 
w. The simulation shown in Fig. [11] shows that even at finite S (non-symmetric initial 
conditions) the final collapse is described by a symmetric solution. 

The second kind is the exploding singularity [153j . since now the two self-similar 
solutions are moving apart, cf. FigJT^ This time even a very small asymmetry 
(5 = 1/100) makes one pinching event "win" over the other. However, this does not 
affect the asymptotics described briefly below. Locally, the solution can be written 

h{x, t) = 6\t')H{{x' - at')/6{t')), (7.7) 

with 6 = t'/ln{t'), which is similar to examples considered in section [31 However, an 
additional complication consists in the fact that the singularity is moving, so there is 
a coupling to the parabolic region between the two pinch-points. This matching is 
unaffected by the fact that in the simulation shown in Fig. [12] one side of the solution 
touches down first. In [153] . a possible generalisation is also conjectured, which has the 
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form 

h{x,t) = S\t')H{{x' - at"^)/6{t')), (7.8) 

r — 1 

and S = t'^~/lnt'. In principle, any value of r is possible, but numerical evidence has 
been found for r ^ 3 (above) and r ^ 5/2 only. 

Finally, a third type is the symmetric singularity of |153j . which does not move. In 
that case, the structure of the solution is 

h{x,t) = hQ{t')H{{x'/6{t')), (7.9) 

with ho = 6^P{ln6), where P is a polynomial. The time dependence of 6 is not reported. 
Evidently, many aspects of the exploding and of the symmetric singularity remain to be 
confirmed and/or to be worked out in more detail. 

The most intriguing feature of the Hele-Shaw equation (17. ip is that several types of 
stable singularities have been observed for the same equation. Within a one-parameter 
family of smooth initial conditions, all three types of singularities can be realized as 
/i — > 0. Each type is observed over an interval of the parameter w. Near the boundary 
of the intervals, a very interesting crossover phenomenon occurs: the solution is seen to 
follow one type of singularity at first (the exploding singularity, say), and then crosses 
over to a solution of another singularity (the imploding singularity). The dynamics of 
each singularity can be followed numerically over many decades in t'. By tuning w, the 
crossover can be made to occur at arbitrarily small values of h. 

The switch in behaviour is driven by the slow dynamics of scaling regions exterior 
to (17. 3p or (17. 7p . It is a signature of the very long-ranged interactions (both in real 
space as well as in scale), that exist in (17. ip . Thus an outside development can trigger a 
change of behaviour that is taking place on the local scale of the singularity. To mention 
another example, applying different boundary conditions for the pressure at the outside 
of the cell can change the singular behaviour completely [154] . This makes the crossover 
behaviour of (17. ip very different from that observed for drop pinch-off (cf. (I2.10p . (l2.1ip ). 
which is driven by a change in the dominant balance between different terms in (12. lip . 

7.2. Semilinear wave equation 

It appears that the Hele-Shaw equation is not an isolated example, but rather is 
representative of a more general phenomenon. Namely, another example of a potentially 
complex singularity structure is the semilinear wave equation 

Utt-Au = \u\P'\ p>l. (7.10) 

It has trivial singular solutions of the form 

u{x,t) = boiT-t)-l^^, (7.11) 

1 

with bo = ^z^Yp- ^ ■ Nevertheless, the existence of different self-similar solutions is 
known in a few particular cases, like the case p > 7, where p is an odd integer (see [155] ) 
or in space dimension d = 1 (see [156] ). 
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The character of the blow-up is controlled by the blow-up curve T{x), which is the 
locus where the equation first blows up at a given point in space. It has been shown 
for d = 1 |157j that there exists a set of characteristic points, where the blow-up curve 
locally coincides with the characteristics of fl7.10p . The set of non-characteristic points 
Jo is open, and T is on Iq. Recently, it has been shown |158] that the blow-up 
at characteristic points is of type II. Even more intriguingly, it appears [158] that the 
structure of blow-up at these points is such that the singularity results from the collision 
of two peaks at the blow-up point, very similar to the observation shown in Fig. [HI 

7.3. More complicated sets 

In the Hele-Shaw equation of the previous subsection, different parts of the solution, 
characterised by different scaling laws, interacted with each other. In the generic case, 
however, finally blow-up only occurred at a single point in space. An example where 
singularities may even occur on sets of finite measure is given by reaction-diffusion 
equations of the family 

Ut-Au = uP-b \Vuf for x e Q. (7.12) 

where Q is any bounded, open set in dimension d. Depending on the values of p > 1 
and g > 1 singularities of fl7.12p may be regional {u blows up in subsets of Q of finite 
measure), or even global (the solution blows-up in the whole domain); see for instance 
[159] and references therein. 

Singularities may even happen in sets of fractional Hausdorff dimension, i.e., 
fractals. This is the case of the inviscid one-dimensional system for jet breakup (cf. 
[160] ) and might be case of the Navier-Stokes system in three dimensions, where the 
dimension of the singular set at the time of first blow-up is at most 1 (cf. [161] ). 
This connects to the second issue we did not address here. It is the nature of the 
singular sets both in space and time, i.e. including possible continuation of solutions 
after the singularity. In some instances, existence of global in time (for all < t < oo) 
solutions to nonlinear problems can be established in a weak sense. For example, this 
has been achieved for systems like the Navier Stokes equations [162] . reaction-diffusion 
equations [163] . and hyperbolic systems of conservation laws [96]. Weak solutions allow 
for singularities to develop both in space and time. In the case of the three-dimensional 
Navier-Stokes system, the impossibility of singularities "moving" in time, that is of 
curves x = if{t) within the singular set is well-known [161] . Hence, provided certain 
kinds of singularities do not persist in time, the question is how to continue the solutions 
after a singularity has developed. 
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